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Abstract : The behaviour of the photon renormalization function in strong 
coupling QED has been recently studied by Kondo, Mino and Nakatani. We 
find that the sharp decrease in its behaviour at intermediate photon mo- 
menta is an artefact of the method used to remove the quadratic divergence 
in the vacuum polarization. We discuss how this can be avoided in numerical 
studies of the Schwinger-Dyson equations. 



As part of a longer study of chiral symmetry breaking in strong QED with Nf flavours we 
have turned our attention to the results of Q where a solution of the simultaneous Schwinger- 
Dyson equations in strong coupling QED is presented in a self-consistent way. Previous studies 
have most commonly approximated the photon propagator by its one loop perturbative form 
in undertaking either analytic § or numerical calculations In contrast, Kondo et al. 
have studied a fully coupled system of equations for the photon and fermion propagators. Then 
the photon renormalization function is determined in a way that is claimed to be self-consistent. 

The Schwinger-Dyson equations for the fermion propagator and for the photon propagator in 
QED are given diagrammatically in Fig. |^. 
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Figure 1: Schwinger-Dyson equations for the fermion and photon propagator. 
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Substituting iSp for the fermion propagator, iD^jy for the photon propagator and (—ieT^) for 
the vertex yields : 
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where q = k — p, and 
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d^k{iT^{k,p))iSF{m^'')iD^,{q) , 



where p = k — q. 

We define the full fermion propagator of momentum p by : 
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The bare fermion propagator is : 
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The full photon propagator of momentum q is given by : 



iD^v{q) = — 2 
The bare photon propagator is : 



0{q ) Qf^u 
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From Eq. (|) one can project out the integral equations for S(p^) and J-{p^). In Minkowski 
space these are given by : 
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Tr[r^(A:,p)(K + S(fc2))^-] 
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Tr[^r^(fc,p)(K + S(A:2)h^ 



It is important to note that unless the vertex T^[k^p) satisfies the Ward-Takahashi identity and 
the regularization of the loop integrals is translation invariant, the photon propagator of Eq. (Q) 
will not have the Lorentz structure of Eq. with the coefficients of g^'^ and q^q'^ being related 
to a single function Q{q^). When these conditions are satisfied then the integral equation for 
Q {q^) can be deduced by applying the projection operator V^u = g^u — nq^qy/q^ (with any value 
of n) to Eq. (pi) : 



1 



a((?2 



iNfe^ 1 
3(2^^ 



d^k 



1 



(/fc2 _ S2(A;2))(p2 _ S2(p2)) 

X P^,TV[r^(A;,p)(K + nk'Wii + S(p2))] . 



(7) 



In general, if we regularize the theory using an ultraviolet cutoff the vacuum polarization integral 
in Eq. ^ contains a quadratic divergence which has to be removed, since such a photon mass 
term is not allowed in more than 2 dimensions. One can show that the q^qu/o' term of the 
transverse part cannot receive any quadratic divergent contribution. Consequently, if we choose 
the projection operator P^jy of Eq. (^) with n = 4, the resulting integral will be free of quadratic 
divergences because the contraction V^^g^^ vanishes. 

A much used alternative procedure is to take the projection operator in its simplest form, 
'Pfj.u = g^lu■ The resulting vacuum polarization integral then contains a quadratic divergence 
which can be removed explicitly by imposing : 

I™ r7^ = 0' 
g2_o g[q^) 
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to ensure a massless photon. If we write the photon renormahzation function as : 

Eq. then corresponds to a renormahzation of the vacuum polarization 11 (g^) : 

q'U{q^) = q^U{q^) - lim q^U{q^) . (9) 

This is the procedure adopted by Kondo et al. [Q]. They solve numerically the coupled set of 
integral equations for the dynamical fermion mass S(p^) and the photon renormahzation function 
in the case of zero bare mass, mo = 0. The calculations are performed in the Landau 
gauge with the bare vertex approximation, i.e. V^^kjp) = As a further approximation they 
decouple the ^-equation by putting T{p^) = 1. While the quadratic divergence in the vacuum 
polarization is removed by imposing Eq. (^), the fact that the Ward-Takahashi identity is not 
satisfied, when dynamical mass is generated, makes the results procedure dependent. 

The integral equations one obtains using these approximations, transformed to Euclidean space, 
changing to spherical coordinates and introducing an ultraviolet cutoff on the radial integrals, 
are given by : 

where (f' = + k"^ — 2pk cos 9, with p = k = \/A?, and 
1 ANfa 1 f^^ „2 A;2 



dk^-, — -TTTiTT / d9 sin^e (111 







g(g2) 37r2 q2Jo A;2 + S2(A:2) 

/ fc^ - A:gcosg + 2S(A:2)S(p2) ^ fc2 + 2S2(fc2) 

^ \ p2 + 5]2(p2) fc2 + 5]2(^2) 

where p^ = q^ + k'^ — 2qk cos 9, with q = k = \fW' . 



The second term in {• • •} in Eq. ( |llD subtracts the quadratic divergence. Recall that in QED the 
momentum dependence of the coupling comes wholly from the photon renormahzation function, 
so solutions for Q{q'^) give the running of the coupling. Kondo et al. solve this coupled set of 
non- linear integral equations, Eqs. (|lO[ 11), for Nf = 1 and find a symmetry breaking phase for 



a greater than some critical coupling ac ~ 2.084. 

In Figs. ^ |3| We display the results for a value of a = 2.086, close to its critical value. The 
dynamical mass function, T,{p'^), is illustrated in Fig. |^. Fig. ^ shows the photon renormahzation 
function, G{q'^), found from their self-consistent solution and this is compared with its 1-loop 
approximation. One observes that at high momenta the self-consistent Q{q^) follows the 1-loop 
result very nicely. For decreasing momenta the effect of the dynamically generated mass comes 
into play and the value of Q{(f'), and hence that of the running coupling, seems to stabilize for a 
while, as one could expect. Then, surprisingly, at some lower momentum there is a sudden fall 
in Q{q^), which drops below the 1-loop value and almost vanishes completely. This is a rather 
strange behaviour for the running coupling at low momenta. This decrease corresponds to tl{q^) 
of Eq. (^) becoming large. 
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Figure 2: Dynamical mass function as a function of momentum for Nf = 1 and 

a = 2.086 as calculated in a self-consistent way as in |||] (A = 10^). 
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Figure 3: Photon renormalization function Q{q'^), as a function of momentum for Nj = 1 and 
a = 2.086 as calculated in a self-consistent way as in [|| and in 1-loop approximation (A = 10^). 
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To solve the problem numerically Kondo et al. have made supplementary assumptions about 
the ultraviolet behaviour of S(p^) and G{q^)- These arise from the need to handle loop momenta 
beyond the ultraviolet (UV) cutoff. For example, if in Eq. ( [To| ) < p^,k^ < A^, then the photon 
momentum q'^ = + k"^ — 2pk cos 9 will lie in the interval < < 4A^ . The same argument 



holds for the fermion momentum p in Eq. dnD, i.e. < < 4A\ As a consequence the 



angular integrals need values of S and G at momenta above the UV-cutoff, this is outside the 
physical momentum region. Therefore one will have to extrapolate S and Q outside this region. 
In their work, Kondo et al. define : 

> a2) = (12) 
n(g2 > a2) = ^ Giq'^ >A^) = 1. (13) 

Both dynamical mass and vacuum polarization vanish above the UV-cutoff and the theory then 
behaves as a free theory. Although this assumption seems reasonable, Eq. ( [T^ ) introduces a 
jump discontinuity in the dynamical mass function at = A^ because ^(A^) 7^ for a > Oc 
(see Fig. |2|), while Eq. (|l^) introduces a relatively sharp kink in the photon renormalization 
function at that point (see Fig. |3|). 

In the physical world these functions have to be smooth. To investigate in a crude way the 
influence of the discontinuity in S(p^), we can remove it by hand by defining the following 
simple extrapolation rule : 

S(/ > a2) = s(A2) ^ . (14) 

This will get rid of the jump discontinuity in the dynamical mass function, leaving instead a 
very slight kink. 

When solving the integral equations using this extrapolation rule, the step in the photon renor- 
malization function at intermediate low momenta surprisingly disappears as can be seen in Fig. ^ 
This was not anticipated since one would not expect the high momentum behaviour of ^{p'^), 
where its value is quite small, to play such a major role in the behaviour of Q{q^) at low q^. 

A more detailed investigation indeed shows that the step in the photon renormalization function 
found by Kondo et al. is an artefact of the way they renormalize the quadratic divergence in the 
vacuum polarization integral, Eq. (pT[), combined with the presence of the jump discontinuity 
in the dynamical mass function, Eq. (^), as we now explain. 



From the angular integrand of the ^-equation, Eq. (11) , we define fe as : 



_ A;2 - fcgcosg + 2S(A;2)S(p2) ^ k'^ + 2T?{k'^) 
■'^ ~ p2 + 5]2(^2) A;2 + S2(A:2) ' ^ ' 

Both terms in Eq. ( [T5| ) cancel exactly at = to remove the quadratic singularity. It is easy 
to see that provided S(A:^) is continuous for all k'^, fg will be continuous, and if S(fc^) has a 
Taylor series, fg will be smooth. Of course the description of the real world has to be such that 
the approximate cancellation of the quadratic divergence at low q^ becomes exact at = in 
a smooth way. 
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Figure 4: Photon renormalization function Q((^), as a function of momentum for N ^ = 1 and 
a = 2.086 as calculated in a self-consistent way with a continuous extrapolation for with 
the jump discontinuity in S(p^) as in |j|] and in 1-loop approximation (A = 10^). 



Now let us look at the angular integrand fg in the approximation of Kondo et al. [|| when 
q'^ is small but k'^ is very large, indeed larger than /cq = (A — g)^. For values of 6 greater 
than 0o(^^) = arccos((A;^ + q'^ — h?')/2kq) we will have > A^. If we now use Kondo et al's 
extrapolation, Eq. ([l^) , then S(p^ > A^) = and the angular integrand Eq. ([Tsl) , now becomes : 

_ k^-kq cos 9 k'^ + 2Y?{k'^) 
is - J2 A;2 + S2(A:2) ' 



When g — > 0, i.e. p ^ k 



fe - -T2f&^+C^(9^^'/cos^). (17) 



As soon as (f' deviates from zero, the angular integrand contains a jump discontinuity at ^ = 
^o(^^)) and part of the angular integrand will not vanish continuously when q^ 0. In fact the 
angular integral Iq will receive an extra contribution 5Io when k"^ is larger than /cq = (A — g)^ : 



dO sin^ e 



T?{k^ 
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k2 + S2(A;2) 
sin 26^0 (A:^)\ 



S2(fc2 



k2 + S2(A;2) 



(18) 



Substituting Eq. (|T8|) in Eq. (y) we see that the vacuum polarization receives an extra contri- 
bution 6IL{q'^) : 



OT(g2) 
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(A;2 + S2(fc2))2 



go(A;2) ^ sin2go(fc2) 



(19) 
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Writing k = A-\- q cos ip, so that — for << A'^ , we have, using the mean value theorem : 

(20) 



, 2^ 8Nfa A3s2(A2) 

^n(? ) - 3,2 ,(A2 + S2(A2))2X/2 ^^^^"^ 



TT ■i/' ^ sin 2^ 



2 2 4 



so that : 



Because of the 1/q this change in n(g2) would be noticeable at very small values of q^. How- 
ever, this analytic calculation does not explain the sharp decrease of G{q'^) at intermediate low 
momenta we and Kondo et al. Q find — see Fig. ^. 

To understand why this happens we have to consider how the numerical program computes the 
extra contribution Eq. (|l9|) to the vacuum polarization integral. The integrals are approximated 
by a finite sum of integrand values at momenta uniformly spread on a logarithmic scale. For 
small g2 the extra contribution is entirely concentrated at the uppermost momentum region of 
the radial integral (fc^ g [fcgjA^]). There the numerical integration program will have only one 
grid point situated in the interval [fcg, A^] for any realistic grid distribution. This point will lie 
at k'^ = Is? if we use a closed quadrature formula. Therefore the integral will be approximated 
by the value of the integrand at A^ times a weig ht factor W{IS?) = wP? {w is 0(1)) : 



4Ar;« Ty(A2)A2s2(A2) e,{K^) sin2go(A2) ^ 

CiiWj-- 3^2 g2(A2 + S2(A2))2 1^2 2^4^ 



(22) 



For small (f' we have 6Q{h?) ~ 7r/2 and the extra contribution to the vacuum polarization will 
be : 

,n(,^).^5!(^, (.3) 

This will effectively add a huge correction to the vacuum polarization at low q^. This has 
been extensively checked numerically and shown to be completely responsible for the sudden 
decrease in the photon renormalization function Q{(f') at low momenta. To reproduce our 
previous analytic result of Eq. (^) numerically, the integration grid has to be tuned unnaturally 
fine to include more points in the region [fcQ,A2]. Without such tuning one has the result of 
Eq. (p3|). Then q^Jl{q^) does not vanish smoothly as q^ — > 0. Instead, for q^ > 0, q'^Il{q'^) 
NfawT,'^{A'^)/3ir and so as soon as g2 is non-zero the cancellation of the quadratic divergence 
disappears suddenly and not gradually as the physical world requires. 
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How can we avoid this problem? As discussed before one can introduce a smooth decrease of 
S(p^) above the UV-cutoff. This ensures that the cancehation of the quadratic divergence takes 
place smoothly as cf' 0. The results obtained with the approximation of Eq. (14) are shown 
in Fig. ^ and are consistent with our physical intuition about the behaviour of the running of 
the coupling. 

Once the quadratic divergence has been removed properly, other numerical difficulties start to 
show up. For instance, inadequate interpolation may give rise to unphysical singularities in 
Q{q^)- We do not discuss these further, as they are outside the scope of this note. However, we 
remark that these problems are avoided if one uses some smooth solution method. 

We conclude that one has to ensure the proper removal of the quadratic divergence from the 
vacuum polarization integral when solving the coupled set of integral equations for the dynami- 
cal mass function and the photon renormalization function numerically. As shown, a very small 
jump discontinuity in the extrapolation of the dynamical mass function can alter the behaviour 
of the photon renormalization function quite dramatically at low momentum and such a peculiar 
running of the coupling is unphysical. To avoid this and also other numerical problems encoun- 
tered in the solution of the coupled set of integral equations it would therefore be preferable to 
search for smooth solutions for the dynamical mass function the fermion wavefunction 

renormalization J-{p'^) and the photon renormalization function G{q'^)- A study implementing 
this is currently in progress. This is essential if we are to understand the phase structure of 
strong coupling QED in 4 dimensions in the continuum. 
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